Analysis of atomic-orbital basis sets from the projection of plane-wave results 
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The projection of the eigenhinctions obtained in standard plane-wave first-principle calculations is 
used for analyzing atomic-orbital basis sets. The "spilling" defining the error in such a projection al- 
lows the evaluation of the quality of an atomic-orbital basis set for a given system and its systematic 
variational optimization. The same projection allows obtaining the band structure and the Ifamil- 
tonian matrix elements in the previously optimized atomic basis. The spilling is shown to correlate 
with the mean square error in the energy bands, indicating that the basis optimization via spilling 
minimization is qualitatively equivalent to the energy-minimization scheme, but involves a much 
smaller computational effort. The spilling minimization also allows the optimization of bases for 
uses different than total-energy calculations, like the description of the band gap in semiconductors. 
The method is applied to the characterization of finite-range pseudo-atomic orbitals [O. F. Sankey 
and D. J. Niklewski, Phys. Rev. B 40, 3979 (1989)] in comparison to infinite-range pseudo-atomic 
and Slater-type orbitals. The bases are evaluated and optimized for several zincblende semiconduc- 
tors and for aluminum. The quality of the finite-range orbitals is found to be perfectly comparable 
to the others with the advantage of a limited range of interacions. A simple scheme is proposed 
to expand the basis without increasing the range. It is found that a double-z basis substantially 
improves the basis performance on diamond, whereas d-polarization is required for Si and Al for 
similar results. Finally, the projection allows the chemical analysis of the plane-wave results via 
population analysis on the previously optimized atomic basis. The spilling can partially control the 
intrinsic arbitraryness of the population analysis. 



I. INTRODUCTION 



One of the several approximations needed for electronic 
structure calculations of solids is the truncation of the 
one-particle Hilbert space. Choosing an appropiate basis 
is critical for obtaining high-quality results. There are 
schemes based on localized, extended, or mixed bases. 
Most of the methods of the first group use atomic-like 
bases, the eigenfunctions being obtained as liaear combi- 
nations of atomic orbitals (LCAO methods) JJl The mosat. 
widely used among the extended is the plane- wave basis.B 

Plane waves (PW) together with pseudopotentials 
have shown to offer a very successful calculation scheme 
for a very large number of applications. Particpjlary since 
the development of ah initio pseudopotentialscl PW have 
provided the tool foLvery accurate calculation of different 
properties of solids a PW are quite simple to implement 
and the convergence of the calculations can be controlled 
with a single parameter, the planewave cut-off. The PW 
basis is also independent of atomic positions, which is 
very convenient for the coding. However, in addition 
to these advantages, plane waves have important draw- 
backs, namely (i) the imposed translational invariance 
(supercells for non-periodic systems), and (ii) their inef- 
ficiency in basis size. 

LCAO methods are much more efficient in the size of 
the required basis. This is a very important advantage 
for the calculation of large systems. Moreover ,|-they have 
shown to be very suitable for order-N methods! in which 
the computational effort scales linearly with the size of 



the system. However, this large reduction in the size of 
the basis is obviously accompanied by a potential loss of 
completeness which can affect the results. It is then nec- 
essary to choose an appropiate basis to obtain accurate 
results. Nevertheless, it is important to stress that, once 
a high-quality atomic basis set has been chosen, calcula- 
tions can be performed with an accuracy perfectly com- 
parable to that of plane waves and with a large reduction 
of CPU time and memory. 

Several methods have been developed to optimize 
LCAO basis sets. Most of them are based in two pro- 
cedures: (i) mininiization of total energies of atoms, 
molecules, or solids,El and [ii) minimization of differences 
in the energy bands cumparing with the experiment or 
with PW calculations oEI We present here a projection 
method that links PW and LCAO schemes and pro- 
vides: (i) a systematic way to characterize the quality 
of an atomic basis and therefore a way to optimize it, 
{ii) LCAO energy bands and Hamiltonian and overlap 
matrix elements, and {Hi) chemical information through 
population analysis. The key of the optimization method 
is to measure the ability of a set of atomic orbitals to 
represent the actual eigenstates of the PW calculation 
of a system. By minimizing the error in the projection 
of the eigenstates the atomic basis can be optimized for 
that particular system. This scheme is related to meth- 
ods of quantum chemistry where a relatively small LCAO 
basis set is optimized by minimizing the distance (max- 
imizing the overlap) of its subsj>ajse to the one spanned 
by another, larger LCAO basisflo The scheme used in 
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this paper was already outlined in Ref. ^, the projection 
technique for bands being already sketched by Chadi in 
Ref. |. 

Besides the projection optimization, the possibility of 
calculating LCAO band structures and matrix elements 
allows the characterization of (i) the energetics, and (ii) 
the range of interactions of the basis. Particulary for this 
purpose non-orthogonalized atomic basis sets are more 
adequate, since the interaction parameters have proved 
to be much more transferahlp and to have shorter range 
than for orthogonal bases.EJO The direct use of atomic 
orbitals in condensed systems gives non-orthogonal basis 
sets naturally. These non-orthogoiaal atomic basis sets 
will be used throughout this paper .E3 

In the present paper, this method of analysis and 
optimization is applied to the characterization of ba- 
sis sets made of the finte-range pseuda-atomic orbitals 
introduced by Sankey and NiklewskiJlj They are coifc 
pared with the infinite-range pseudo-atomic orbitalstil 
and Slater- type orbitals. Both minimal and expanded 
bases are explored as applied to several zincblende semi- 
conductors, as well as to aluminum. 

The optimized basis can be used for the calculation 
of larger systems, where plane- wave calculations are im- 
practical. Besides that possible use, they can also be used 
directly for the analysis of the PW results themselves.El 
The chemical language of population analysis is not ac- 
cesible to the PW scheme. The projection into the op- 
timum LCAO basis allows its application, to compute 
atomic charges, bond charges, and charge transfers. In a 
similar way as the decomposition of the total charge into 
orbital populations, the density of states can also be pro- 
jected into atomic orbitals to give local densities of states. 
It is shown here how the optimization of the basis is im- 
portant for obtaining sensible values of the charges, and 
also how the projection can be used to partially control 
the intrinsic arbitraryness of any population analysis. 

The structure of this paper is as follows. A brief de- 
scription of the method to characterize and optimize the 
atomic basis is given in section II, as well as a descrip- 
tion of the way of performing the population analisys 
and the band structure calculations. Section III describes 
the main characteristics of the ab initio PW calculations 
used here. A description of the different types of atomic 
orbitals analyzed in this paper can be found in section 
IV. Section V proposes the spilling as a parameter to 
characterize basis sets and their correlations with other 
quantities. Section VI contains the results of the analisys 
and optimization of the differents sets of atomic orbitals. 
Population analysis is discussed in section VII. We end 
with the conclusions in section VIII. 



II. DESCRIPTION OF THE METHOD 

The method pr^epted in this paper is based on pro- 
jection techniques .Bo Given the results of a PW calcula- 



tion, the projection of the calculated eigenstates into an 
atomic basis can provide information about this basis set 
and about the PW results. A good basis must be able 
to represent the essential features of the PW eigenstates. 
Once the PW calculation has been performed for the sys- 
tem of interest, the projection process is much less costly 
than the self-consistent calculation itself, and many dif- 
ferent basis sets can be analyzed and optimized for one 
single PW calculation, a much more economic procedure 
than performing a full self-consistent calculation for each 
trial basis set. 

On the base of our method we have the definition of the 
spilling parameter.EI Given a PW calculation we define 
the spilling Sn for a given local basis set as 

^^^ = w,n:^ E(^"(k)i(i - p(k))i^a(k)), (1) 

k a — 1 

where |'0Q(k)) are the PW calculated eigenstates, and 
Nk and are the number of calculated k points in 
the Brillouin zone and the number of bands considered, 
respectively, and n is the specific subspace spanned by 
the eigenstates included in the sum (defined by Nk and 
Na). -P(k) is the projector operator into the subspace of 
Bloch functions of wave vector k, generated by the atcffliic 
basis, and defined as usual for non-orthogonal basis,t3 

p(k) = I'p.mirm - \M^))s;,\k){M^)i 

(2) 

where 

(r|,/>^(k)) = ^ <^^(r - - R)e*(---+^) (3) 

R 

(j)^{r) being the atomic orbitals, the atomic coordi- 
nates in the unit cell, and R the lattice vectors. 

^^.(k) = {(bMM^)) (4) 

is the overlap matrix of the atomic basis. The vectors 
|(/)^(k)) constitute the dual of the atomic basis, that sat- 
isfy 

irikM.m = = S^,. (5) 

S(i is a parameter that varies between and 1 and mea- 
sures the ability of a basis to represent PW eigenstates, 
by measuring how much of the subspace of the Hamilto- 
nian eigenstates falls outside the subspace spanned by 
the atomic basis. If we consider the PW eigenfunc- 
tions iV'a(k)) and their projection into the atomic ba- 
sis P(k)|'0Q,(k)), then Sn gives the average of ||(1 — 
P)|V'Q(k))|p over the eigenstates considered for the pro- 
jection. It must be also pointed out that the projected 
eigenstates |xQ,(k)) = P(k)|'(/;Q(k)) are not necessarily or- 
thonormal, the overlap matrix i?a^(k) — (xQ(k)[x/3(k)) 
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being different from the identity matrix. In particular 
RaaO^) < 1, and 

" k Q=l 

is the averaged norm which has been lost (spilled) in the 
projection process. If we restrict our analisys to the occu- 
pied eigenstates, then Sq is the fraction of total electronic 
charge spilled when projecting, and will be refered to as 
charge spilling Sq . 

A local spilling can also be defined as a function of real 
space coordinates, 

Snir) = ^^EE(l'^"(k,r)p - |x.(k,r)n (7) 

" k a=l 

We consider it a useful and direct way for visually charac- 
terizing the quality of the representation of the P W eigen- 
states. The total spilling defined above is recovered by 
integrating the local spilling over the three-dimensional 
space. 

Important for the characterization of the quality of a 
basis for a given system is the energy band structure. 
To calculate the LCAO energy bands of a system we 
project the Hamiltonian associated to the self-consistent 
PW charge density of our reference calculation into the 
basis subspace: 

Hi:^^0{\.) = (0^(k)|ff^>^i0.(k)) (8) 

To obtain the LCAO Hamiltonian matrix elements we 
expand the Bloch basis functions into plane wavesfl We 
calculate the kinetic energy and the effect of the non-local 
part of the pseudopotential directly in Fourier space. The 
effect of the local part of the potential and self-consistent 
Hartree and LDA interactions is calculated using fast 
Fourier transform algorithms, 

^^.^^°(k) = E (^^(k)|k + G)(k + G|ff^^|0.(k)) 

|k-i-G12<£;„„^ 

(9) 

5^.^^°(k)- E (0^(k)|k + G}(k + G|^.(k)) 

|k+G|2<£;„„, 

(10) 

where G are reciprocal lattice vectors and Emax is an en- 
ergy cut-off independent of the one used in the PW calcu- 
lation, that must be large enough to guarantee a reliable 
representation of the atomic basis. With this method we 
obtain a LCAO Hamiltonian from a first-principles PW 
calculation, no free parameters being fitted. 

The LCAO Hamiltonian is obtained in Bloch space 
and, as a consequence, it includes interactions up to infi- 
nite neighbors. The deviations in the band structure with 



respect to the PW only reflect the uncompleteness and in- 
adequacy of the basis, clearly separated from the possible 
additional approximation of neglecting matrix elements 
beyond some range of interactions. We can study both 
approximations separately. In most studies about LCAO 
band structures these two effects are mixed together. As. 
an exception it must bepCnentioned the works of ChadicI 
and Jansen and SankeyJiJ where also LCAO Hamiltoni- 
ans with interactions up to an infinite range of neighbors 
are obtained. 

The Hamiltonian matrix in real space is calculated 
from the Hamiltonian matrix in Bloch space by perform- 
ing a Bloch type transformation 

if,^,^^0(R,.) = E i/,^,^^0(k)e*(^-^'^' (11) 

k 

S^^^^iR,.) = E 5^.^^°(k)e^''(^-^'^) (12) 

k 

where normalization factors which depend on the over- 
laps are omitted for clarity. The sum has to be extended 
to a sufficient number of k points, taking into account 
that the number ofpcitDts depends on the real space range 
of the interactions.tjllj 

The projection technique can also be used for obtaining 
chemical information from the PW calculations by means 
of an LCAO population anahsys. We use the one pro- 
posed by Mulliken.ll3 The PW occupied eigenstates are 
projected onto the subspace spanned by the atomic ba- 
sis. Due to the non-orthogonality of the projected eigen- 
states, the density operator has to be properly defined to 
ensure the charge conservation 

occ 

P = EElXo(k))(x"(k)|, (13) 

k a 

where |x"(k)) — X]/3 ^/3cl(k)|X/3(k)) represent the vec- 
tors of the dual set of the projected eigenstates, and 
i?/3a(k) — (x/3(k)|xQ(k)). The LCAO density matrix is 
then written in terms of the dual LCAO basis: 

r^. = {r\p\r), (14) 

the charge associated to an orbital /i being 

Qf-E^A'-V- (15) 



III. REFERENCE CALCULATIONS 

The analysis of the bases presented in this paper are 
based on reference PW calculations which have been 
performed using ab initio pseudopotentiais generated 
by the improved TrouUier-Martins methodlla and within 
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the local density approximation (LDA) for electron ex- 
change and correlation.Ej Of course, other schemes could 
have been equally used, the only requirement to apply 
the method being the planewave expansion of the one- 
electron wavefunctions. The energy cut-off for the plane- 
wave expansion has been taken large enough to consider 
that all the calculations are well converged. High en- 
ergy cut-offs are not necessary for the application of the 
method but are important to ensure accurate results for 
the optimization of the orbitals as well as for the pop- 
ulation analisys. The cutoffs used are: 10 Ry for Al, 
20 Ry for Si, BP, and AlP, 30 Ry for GaAs, 40 Ry for 
graphite, 50 Ry for BN, and 70 Ry for C and SiC. For the 
self-consistent zincblende semiconductor calculations, we 
used 2 irreducible-Brillouin-zonc points (IBZP), equiva- 
lent to 32 whole-Brillouin-zone points (BZP), giving a 
k-mesh length cutofjtB Ic ~ a, where a is the lattice con- 
stant. For aluminum and graphite we used respectively 
(10 IBZP, 125 BZP, ^^==13.435 a.u.) and (8 IBZP, 16 
BZP, ^£=4.648 a.u.). The length cutoff Ic is simply re- 
lated with the maximum range of the real-space matrix 
elements obtainable from the PW calculation (Eqs. 11 
and 12) by |R^,y| < Ic- For the pseudopotentials, in all 
the calculations we have used the separable, fully pon- 
local formulation following Kleinman and BylanderEj 



IV. ATOMIC BASIS SETS 

A brief description follows of the sets of atomic orbitals 
which are used throughout this work. The studied basis 
sets can be classified in two groups: (i) those which are 
provided by the numerical solution of the atomic prob- 
lem (with the same atomic pseudopotential as used in 
the PW calculaiiniij) , hereafter called PAO's for pseudct 
atomic orbitals,yMand (ii) Slater-type orbitals (STO).eI 

Within the first group, PAOoo are the atomic solutions 
under real atomic boundary conditions (at r = oo). This 
family of bases allows for a scaling factor for the radial 
part acting as a variational parameter to optimize the 
atomic orbitals via the spilling minimization, 

0r(r) = Af0f"(A,r). (16) 

PAOoo represents quite a natural choice, since access 
to PW calculations is usually accompanied by access 
to atomic calculations under the same approximations. 
They provide a good description of the charge density, 
and have already been used in ab initio calculationstfl 
giving accurate results. For covalent materials, with 
small charge transfers, they provide a very good de- 
scription of the pseudo-core region. This is worth con- 
sidering since the PW eigenstates are expanded in the 
atomic basis maximizing the projection without mak- 
ing any distinction between the pseudo-core region (zone 
near the atomic nucleus) and the valence region (inter- 
atomic zone). 



The main disadvantage of the PAOoo 's is the very long 
range of the interactions and overlaps they originate. 
Getting a shorter and perfectly defined range of inter- 
actions is the reason to work with PAO^^. They are the 
solutions of the atomic pseudopotential when forced to 
vanish outside a cutoff radius rc, keeping-jthe continuity, 
as introduced by Sankey and Niklewskilj (they will be 
denoted by specifying the value of rc in atomic units). 
The cutoff radius represents a variational parameter for 
this type of bases. These orbitals seem promising candi- 
dates in calculation techniques for very large systems .□ 

Within the Slater orbitals, we call STOl the con- 
ventional orbitals r"~^ exp"'^'', where n is the principal 
quantum number and the scale factor (3 is taken to be 
the same for all the atomic-orbitals within the same shell; 
ST02 will be the same, but allowing for a different /3 for 
each orbital, (3s Pp Pd\ and ST03 are the orbitals 
r"exp~'^'', where a becomes another variational param- 
eter, no longer restricted to be an integer, and allowing 
both a and /3 to be different for each orbital. 

STO's have been extensively used in LCAO calcula- 
tions in quanturn chemistry and have proven to provide 
accurate results.S They present the additional advantage 
of their very simple mathematical properties. To sim- 
plify the computational task, it is a common practice to 
expand the STO's orbitals as linear combinations of gaus- 
sian functions. These are the so-called STO-NGB orbitals 
where N stands for the number of gaussians orbitals used 
per STO. We have done some tests to compare the results 
obtained using STO and ST0-4G orbitals. No apprecia- 
ble differences have been found in the results. Spilling 
values and optimum exponents arc essentially the same 
in both cases, and no further investigation has been pur- 
sued concerning STO-NG. 

STO orbitals do not seem a priori to be very well suited 
for the description of the pseudo-core region for the pro- 
jection, possibly leading to worsen the spilling at the va- 
lence region while trying to adapt the pseudo-core. To 
explore that possibility we consider a last kind of orbital 
functions, the pseudo-STO, which represent a link be- 
tween the STO and the PAO families. A pseudo-STO 
is a Slater-type orbital outside a certain radius, and the 
numerical solution of the pseudopotential inside that ra- 
dius, continuity of the function and its derivative being 
imposed. 



V. SPILLING ANALYSIS 

This section is devoted to show the adequacy of an 
analysis based on the projection and the spilling for the 
evaluation of a basis set. The total spilling provides a 
natural parameter that can be used straightforwardly 
to measure systematically the quality of a basis. The 
optimization in terms of spilling provides a variational 
scheme, since the minimum attainable is zero, which 
means a perfect basis. 
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Both the projected charge density and the local spilling 
offer a usehil tool for visually characterizing the qual- 
ity of a basis. This is illustrated in Fig. 1, where a 
charge-density contour map for SiC is shown comparing 
the PW results with the results of the projection on the 
PAOoo, showing in addition the local spilling along the 
bond axis. We see that some of the bond structure is lost, 
the charge density being much more spherical around the 
atoms, when projecting on this basis. This can be also 
seen in the fact that the two maxima in the density near 
the carbon nucleus are of approximately the same height, 
while in the PW calculation there is more density at the 
bond side. All these features can be seen in the local 
spilling, which is negative in the regions with an excess 
of projected charge and positive in the zones with a de- 
fect of projected charge. The local charge spilling is neg- 
ative in the proximity of both core regions and positive 
in the bond region. Some of the bond charge goes to an 
anti-bonding region. This feature also tends to make the 
charge-transfer lower for this PAOoc basis (sc!c below). 

For the spilling analysis to be useful, the spilling has 
to be correlated with the quality in the energetics of the 
system. This is tested in Fig. 2. where the mean square 
error of the projected valence bands of crystralline sili- 
con, as compared to the PW ones, is plotted versus the 
charge-spilling Sq, for different families of bases. It is 
calculated without shifting the bands, the absolute zero 
being unambiguously defined (the Hamiltonian operator 
is not changed). This quantity can also be taken as an 
indicator of the error in the total energy due to the ap- 
proximation of the basis. 

The clear correlation between the spilling and the error 
in the energy bands displayed in Fig. 2 indicates that a 
basis optimization based on a minimization of the spilling 
is parallel (at least approximately) to an energy mini- 
mization. This parallelism is even clearer within a basis 
family, where the correlation is much more pronounced 
(different points within the same family are obtained by 
varying the characteristic parameters presented above). 

To finish this section, we introduce a lower bound for 
the spilling. The dependence of Sq on the atomic basis 
enters through the Fourier coefficients of the radial wave- 
functions. Therefore, <Sn can be regarded as a function of 
these Fourier coefficients (/>i,,j(|k -I- G|) which are differ- 
ent for each atom i, orbital /i and each modulus |k + G|. 
Considering all these Fourier coefficients as free varables 
the minimization of <Sn gives a lower bound for Sn- Val- 
ues of this lower bound will be shown as Min. <S in the 
tables. 



VI. ANALYSIS AND OPTIMIZATION OF 
ATOMIC BASIS SETS 

The quality of the LCAO bases presented above is mea- 
sured and optimized for different solids in this section us- 
ing the spilling analysis, and, simultaneously, the LCAO 



band structure is contrasted with the results of plane- 
waves. The results are displayed in Tables I through IV. 
The resulting numbers are abundant since the number 
of possibilities is large. For the different kinds of bases 
introduced in section IV there are different possible basis 
sizes, single-2 or double-z, meaning one or two orbitals 
per atomic symmetry. There is the additional possibility 
of polarizing the basis, i.e., introducing atomic orbitals 
of a higher angular momentum than needed, which, in 
the cases of this paper represents the addition of the d 
orbitals to the minimal sp basis. 

The optimization of the different bases introduces new 
degrees of freedom since their optimization is done with 
respect to different variational parameters (see section 
IV), and also the spilling can be defined for different num- 
bers of bands, depending on the particular application 
the basis is optimized for. 



A. Minimal basis sets 

Table I shows the spillings of the projections onto dif- 
ferent minimal bases for silicon, diamond, and aluminum, 
together with the mean square error in the LCAO bands 
as compared with the PW ones. Only sp bases are con- 
sidered, except for Al, for which also spd bases data are 
shown due to the poorness of minimal sp bases for this 
solid. The numbers displayed for STO bases have been 
obtained after minimization of the charge spilling (the 
spilling of the valence eigenvectors) for the correspond- 
ing parameters. This procedure is illustrated in Fig. 3 
for the specific cases for silicon and carbon, where the 
charge spilling is shown as a function of the STO expo- 
nents. PAO^* stands in Table I for the PAOoo optimized 
via a scale factor. The PAO^^ are optimized with respect 
to the radius of the basis functions,rc, as illustrated in 
Figs. 4, 5, and 6 for silicon, diamond, and aluminum, 
respectively. In Table I, cspccific values of Vc arc chosen. 

All the bases displayed in Table I perform comparably 
well in what concerns the charge spilling. Clearly bet- 
ter does sp for C than for Si and Al, as expected since 
d orbitals play a lesser role for C than for Si and Al. 
The spilling for eigenstates including conduction bands 
is larger than for valence bands only. For sp bases there is 
a factor of 10 to 20. The two main reasons for that are, (i) 
the bases have been optimized by minimizing the charge 
spilling, and (ii) higher energy bands require higher en- 
ergy atomic orbitals. 

Table I also shows the populations of s, p, and d or- 
bitals (if present). It can be seen that the sp^ hybridiza- 
tion is better accomplished in C than in Si (a perfect 
hybridization would show Qg = 1 and Qp = 3). Qd in Al 
is quite appreciable. 

The mean square error in the energy bands is displayed 
as (Ty for the valence bands and for the first two con- 
duction bands ((J1+2 for first two bands, and as+i for 
third and forth, in the case of aluminum). The valence 
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bands are quite well described with sp bases for Si and C, 
but not for Al. The metallicity in Al demands a larger 
contribution of d orbitals, what makes them necessary 
for any quantitative LCAO description of Al. The con- 
duction bands are much worse described by the minimal 
bases shown in the table. 

The dependence in the number of bands to be opti- 
mized through the spilling is shown in Table II for Si, 
C, BN, and Al. Each kind of basis has been optimized 
by minimizing (i) the charge spilling and (ii) a spilling 
including conduction bands (iSg for Si, C, and BN, and 
^4 for Al). It can be seen how an improvement in 53(4) is 
accomplished when optimizing the basis for it, but pay- 
ing the prize in the worsening of Sq. In addition to the 
S and (T values, the band gap region description is eval- 
uated by showing the values of the F-direct and indirect 
gaps, and the position of the minimum of the conduc- 
tion band. An important and systematic improvement is 
observed in these quantities when the optimization is ex- 
tended to unoccupied states. The displayed information 
offers the necessary information for the evaluation of the 
trade-off in the quality of the different spectral regions for 
the choice of a basis for a particular application. Fig. 7 
shows the band structure of Si obtained with a minimal 
sp ST03 basis optimized considering Sg ■ Notice the nice 
agreement in the band-gap region, showing a very realis- 
tic indirect gap, quite rare for a minimal basis [compare 
with the typical cases shown in Fig. 8 (a) and (d)]. The 
LCAO bands have been shifted by 1.04 eV, what means 
an important deviation in total energies, if computed. 
Another price to pay is the very long range interactions 
due to the unusual extension acquired by the orbitals in 
the optimization. 

Trying to improve the performance, other, more com- 
plicated, kinds of bases have also been tried, the results 
not being shown in the Tables. Pseudo-STO's (see sec- 
tion IV) provide good bases, but the improvement on the 
others (STO and PAO) have not been substantial enough 
for seriously considering them as alternative candidates. 
They have been optimized by varying the STO parame- 
ters outside the core region (as for ST03) plus the extra 
degree of freedom given by the core radius. For Si a mini- 
mum in Sq is obtained for a core radius of 1 a.u., but the 
spilling is only 4 % lower than the spilling of the ST03 
case. This fact indicates that the projection is not as 
sensitive to the core regions so as to worsen the overall 
projection if the core region in the basis orbitals is not 
adapted to the pseudopotential. 

Another kind of basis has been constructed based on 
PAOr^ . More degrees of freedom are given to the spilling 
minimization by linear combining excited states of the 
atomic problem to the ground state. For each symmetry 
{s,Px,Py,Pz, ■■■) the atomic orbital is constructed with a 
linear combination of the ground and excited states of 
the corresponding pseudopotential. This procedure con- 
siderably reduces the spilling without increasing the size 
of the basis. However, the band structure features are 
not accordingly improved. This fact leads to the conclu- 



sion that, for the spilling-energy correlation to be satis- 
factory in the solid, the energetics in the basis itself (in 
the atomic problem) has to be considered. 



B. Range of interactions 

An important aspect defining the applicablity of 
LCAO bases is their range of interactions, i.e., the scope 
of neighbor shells that have to be included in the LCAO 
Hamiltonian and overlap matrices. The basis sets giv- 
ing comparable results in the previously shown quantities 
perform quite differently in this respect. This is shown in 
Table III. This information is obtained by computing the 
atomic matrix elements of the Hamiltonian and overlap 
matrices as described in section II. Due to the different 
analytical expressions for the orbitals, to measure the in- 
teraction range we adopt an empirical criterium, namely, 
the radius r99.9% for which 99.9 % of the norm of the 
orbital is within the sphere with that radius. We have 
found that criterium meaningful in the sense that neglect- 
ing interactions beyond that radius does not appreciably 
alter the energy bands corresponding to the valence and 
some conduction bands (average error lower than 0.05 
eV). Neglecting interactions within that radius produces 
appreciable deviations in the bands. For PAO^^'s the in- 
teracting scope is clearly defined by Tc (which incudes 
100 % if the norm within), but for short Tc's they are 
very similar to ^99 9%, and for larger re's then r99 9% be- 
comes more meaningful. Therefore, we also characterize 
the PAOr^'s with ^99 9%. Besides this quantity. Table III 
also shows the scope of non-negligible neighboring shells 
for each case. 

PAOoo represents the best basis in what spilling 
and band quality is concerned, and also for population 
purposed^ (see next section), but its range of interac- 
tions is very large. Xhis fact was already pointed out by 
Jansen and SankeylU in their work doing self-consistent 
calculations using this basis. The results of our calcula- 
tions show that interactions up to at least 5 (10) shells of 
neighbors have to be kept for Si (C) to obtain a reason- 
able band structure with the correct semiconducting be- 
havior. The best basis suited for a shorter range of inter- 
actions together with keeping good standards in spilling 
and band structure is the PAO^^ , as can be seen in Table 
III and Fig. 4 for Si, in Fig. 5 for C, and in Fig. 6 for Al. 
Note that, in spite of C needing a smaller Tc than Si to 
attain small spillings, the range of neighbors is larger for 
C due to the much smaller lattice parameter. 



C. Polarized and double-^ basis sets 

The clearest way of improving the quality of a basis is 
increasing its size. In LCAO there are two usual ways of 
doing it, (i) increasing the number of orbitals within the 
same atomic symmetries as found in the minimal basis. 
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and (ii) polarizing the basis. Within the first, double- 
z is the most usual case, i.e., doubling the basis, which 
in the cases of this paper means a basis of two s ,two 
Pa;, two Pj, and two pz orbitals. A polarization of the 
basis is accomplished when a shell of different atomic 
symmetry is added to the minimal basis. In our cases it 
means the addition of the d shell to the minimal sp basis. 
Combination of these two procedures is also customary 
for further improvement. A rule of thumb in quantum 
chemistry says that a iasis should always be doubled 
before being polarized.B This is not the usual practice 
within the LCAO calculations in the solid-state-physics 
community, where double-z bases are hardly to be found- 
Expanded STO bases can be found in the literature. eI 
This is not the case for PAO's. In this paper we propose 
to use the excited states of the atomic calculations lead- 
ing to the PAOr^'s, under the same boundary conditions, 
to obtain the orbitals required to expand the basis, i.e., 
excited s and p orbitals for double-z bases and d orbitals 
for polarization. The advantages of this scheme are its 
simplicity and the controlled range of interactions also for 
the added orbitals, which is usually lost in the traditional 
scheme. 

Again, we study here the quality of expanded LCAO 
bases by means of their spilling and their band struc- 
ture as obtained from the projection. The results of our 
analysis are shown in Table IV, where minimal (singie-z), 
double-z, and polarized single-z bases are compared for 
silicon and diamond. See also Table I for aluminum. The 
sp and spd bases are optimized as above using ST03 for 
the Slater orbitals, and using the corresponding Tc for the 
PAO's, as shown in Figs. A-^ Double-z for STO has been 
taken from the literatureBONote that, with double-z. Si 
does not reach equivalent quality as C, and that, for that 
purpose, the polarization of the basis is quite enough. 

The band structures corresponding to the data of Table 
IV are shown in Fig. 8 for Si, Fig. 9 for C, and Fig. 10 for 
Al (Table I), in each case comparing the LCAO bands ob- 
tained from the projection with the reference PW bands. 
The information given in the Table is ratified, and ob- 
served in more detail. Note that the substantial improve- 
ment in the Si case occurs when including the d orbitals, 
but not when doubling the minimal basis. For diamond, 
however, such qualitative improvement is given by the 
doble-2; basis. The figures for the gap in the Si and C 
cases are shown in Table V. In the Al case, the polar- 
ization is needed from the start, the sp projection giv- 
ing defective band structures in the neighborhood of the 
Fermi level around the W point. The spd band structure 
is shown in Fig. 10 for PAO^^ for different Tc. The perfor- 
mance of the double and polarized STO bases in the first 
conduction bands is better than that of the correspond- 
ing PAOr-^ bases. This is due to the larger extension of 
the extra orbitals in the STO bases, which is not allowed 
in the PAO^-c's. The latter, however, also represent a 
considerable improvement over the minimal bases, while 
keeping the interaction range. This virtue makes them 
very appealing for accurate LCAO calculations. 



The error in the band structures could be noticeably 
reduced with a rigid shift of the bands, which is clear in 
sight of the Figures. This has not been allowed in order 
to keep the band-structure error as representative of the 
behavior to be expected of total energies. 



D. Geometry and environment dependence 

The projection procedure proposed in this paper is 
based on PW calculations of actual materials. As com- 
pared to basis optimization in atoms, this allows adequat- 
ing the basis to particular environments or geometries. 
As a first example we show the change in the minimal 
sp PAOr^ basis for C depending on the solid being dia- 
mond or graphite. Fig. 11 shows the charge spilling as a 
function of and of the range of interacting neighboring 
shells. Graphite is worse represented by the same kind 
of basis. In addition, graphite requires larger to ap- 
proach minimum spilling. This latter fact is due to the 
comparatively large interplanar distances in the graphite 
structure. 

Another example of geometry dependence of the basis 
quality is given in Fig. 12, where the charge spilling of 
a minimal PAO^^ for Si is plotted versus the lattice con- 
stant for different values of Tc- From the Figure it is clear 
that (i) for spilling purposes and, therefore, for the total 
energy, the basis is better the larger rc, (ii) the slope in 
Sq versus lattice constant will translate in a correspond- 
ing slope of the total energy, what originates a fictitous 
pressure on the system, (iii) this pressure is different de- 
pending on Tc, and the optimum for this purpose (the 
one giving zero pressure, with Tc between 5 and 6 a. u.) 
is different from the optimum for energy considerations 
(tc ~ oo), and (iv) there are no appreciable curvatures, 
from where it can be inferred that no substantial devia- 
tions are to be expected in the ellastic and force constants 
(vibrations) due to the basis. 



VII. POPULATION ANALYSIS 

One of the uses of the projection into LCAO is the 
possibility of performing population analysis of the PW 
results, otherwise inaccesible. It is well known that any 
population analisys has an inherent arbitraryness associ- 
ated to the fact that the relevant quantities in the anal- 
ysis do not correspond to observables, but rather to con- 
venient ways of dissecting the total charge of the system. 
Furthermore, small variations in the basis can produce 
quite different population results, particulary for charge 
transfers. In spite of it, much of the chemical descriptions 
of condensed systems is performed in the language given 
by such population analyses. We, therefore, consider it 
worthwhile to explore the problem with the spilling. 

Due to the mentioned arbitraryness, the correlation 
of charge transfers and charge spilling is expected to be 
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poor. This is clear when uncompensating the basis, Hke, 
for instance, when introducing d-orbitals for sihcon in 
SiC. This greatly improves the description of the eigen- 
functions and reduces the spilling, but being the basis 
uncompensated, some features (structure of the bond) 
are described by the much more complete basis of one of 
the atoms, so that more electronic charge is assigned to 
the atom with more complete basis, getting an unphys- 
ical charge transfer. A compensated basis is therefore 
important for the populations to be meaningful. A mea- 
sure of this compensation can be obtained if the spilling 
is calculated independently for structurally similar solids 
containing the individual atoms separately. The spillings 
for silicon and diamond separately, using the same kind 
of basis for both, indicate a slight uncompensation fa- 
voring diamond, which, in this case says that the charge 
transfer might be slightly overestimated. In addition to 
that, the optimization of the basis is important to get re- 
liable values of charge transfers. This is seen in Fig. 13, 
where the charge transfer in SiC is shown to correlate 
with the spilling. A 'converged' charge transfer is not 
guaranteed, but, provided that the basis is compensated, 
a small spilling indicates a meaningful charge transfer. 

Table VI shows the charge transfers for some 
zincblende semiconductors, obtained with the projec- 
tion procedure, using optimized PAOoo basis sets. The 
numbers compare well with equivalent resultSr-abtained 
from direct LCAO selfconsistent calculations .E3'E3 Sys- 
tems with only one species do not present these difficul- 
ties for the charges since the basis can be trivially com- 
pensated. Table VII shows the excess of charge (with 
respect to the neutral atom) assigned to the upperj-and 
lower atoms of a dimer in the Si(100)-2xl surfaceO A 
minimal basis has been used. The populations are quite 
stable and independent of the basis used. 



the uncompleteness of the basis and to considering a lim- 
ited range of interactions. This is due to the fact that 
the bands are obtained directly in Bloch space, with in- 
teractions to infinite neighbors. 

(iii) Both STO's and PAO's optimized with spilling 
minimization give quite good basis sets, of comparable 
quality. 

(iv) If a controlled range of neighbor interactions is 
desired, the best performing basis sets are pro\rided by 
the PAOr^ introduced by Sankey and Niklewski.ll3 

(v) A simple and interaction-range-controlled scheme 
for extending PAO's basis sets to double-z and/or po- 
larized bases is proposed, and shown to give good basis 
sets. 

{vi) Double-z bases substantially improve the perfor- 
mance on diamond, whereas d polarization is required for 
similar results in Si and Al. 

(vii) The spilling minimization allows optimization 
taylored for particular uses, like for the description of 
excitations or band gaps, environment dependence of the 
basis, etc. 

(viii) Population analysis can be performed, and its in- 
trinsic arbitraryness can be partially controlled by means 
of the spilling. 
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VIII. CONCLUSIONS 

We have presented a projection method which provides 
a tool to characterize and optimize atomic orbital basis 
sets for a given system. Independently of how many bases 
are evaluated for that particular system, only one ref- 
erence self-consistent planewave calculation is required. 
The method also provides LCAO band structures, Hamil- 
tonian matrix elements (non fitted tight-binding param- 
eters), and population analysis. 

Based on that tool and on its characteristic parameter, 
the spilling, an exhaustive analysis of two kinds of LCAO 
basis sets, namely. Slater-type (STO) and pseudoatomic 
(PAO) orbitals, have been presented for several solids, 
arriving to the following conclusions: 

(i) The spilling provides a very convenient tool for 
studying and variationally optimizing basis sets, that 
closely correlates with energy optimizations. 

(a) The LCAO band structures obtained from the pro- 
jection allow to characterize separately the errors due to 
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FIG. 1. Electron charge density for SiC: (a) self-consistent 

density contour map in the (110) plane calculated with PW, 
(b) the same after projecting the cigcnstatcs into a PAOoo sp 
basis, and (c) these two densities (PW with solid lines and 
the projected density with dotted lines) and Sq (r) along the 
bond axis. Units are electrons per unit cell. 

FIG. 2. Mean square deviation in valence bands of the 
projected hamiltonian from the PW bands versus the charge 
spilling of different bases in silicon. Each symbol stands for 

a different family of bases: □ for PAOr^, x for PAOoo, + for 
STO, and O for different double- 2; bases. Bases with ay lower 
than 0.2 eV include d orbitals. 

FIG. 3. Charge spilling versus exponent for a STOl basis 
set, (a) silicon with a spd basis, (b) diamond with a sp basis. 
Also shown the optimum exponents for a ST02 basis set for 
the same systems. 



FIG. 6. The same as Fig. 4 for aluminum. 

FIG. 7. Gap optimization with a minimal basis. Silicon 
band structure calculated with a ST03 sp optimized to re- 
produce the occupied eigenstatcs and the first four conduc- 
tion bands. The PW band structure is shown with dots. The 
ST03 bands have been shifted 1.04 eV to make the top of 
both valence bands coincide. 

FIG. 8. Band structure of Si for minimal, double-^, and 
spd bases. In (a), (b) and (c) the Hamiltonian is projected 
in a ST03 sp, STO double-z, and ST03 spd basis, respec- 
tively. ST03 bases optimized for valence band. (d),(e) and 
(f) correspond to PAO5.0 sp, double- 2;, and spd basis. The 
PW reference band structure is shown with dots. 

FIG. 9. Band structure of diamond for single-2; and 

double-2: bases, (a) and (b) correspond the projection in a 
ST03 sp and STO double-2 basis sets, respectively. ST03 
basis optimized for valence band. Figures (c) and (d) corre- 
spond to PAO4.0 sp and double-2;. The PW reference band 
structure is shown with dots. 

FIG. 10. Aluminum band structure: (a) calculated with 
PAOs.s spd basis, and (b) using PAOe.i. Dotted lines show 
the PW band structure (Fermi energy is 8.2 eV). 

FIG. 11. (a) Charge spilling for diamond (D) and graphite 
(G) as a function of Tc for a PAOr-e *P basis, (b) Charge 
spilling for diamond as a function of the shell of neighbors 
with non-zero overlap, and (c) the same for graphite. 

FIG. 12. (a) Sq versus lattice constant of silicon for PAOr^ 
bases, (b) SSq—Sq-Sq where the Sq is taken at the equilib- 
rium value of the lattice constant oo=10.2 a.u.. O stands for 
rc=4 a.u., + for rc=5 a.u., □ for rc=6 a.u., x for Tc = 00. 

FIG. 13. Charge transfer in SiC as a function of the charge 
spilling of the sp basis used to project the density. □, stands 
for PAOoo optimized for the heteropolar system and indepen- 
dently for Si and C in a zinc-blende structure. O, stands for 
ST03 basis in the same conditions, and -|- for a ST02 basis. 



FIG. 4. (a) Charge spilling for silicon with a PAOrc *P 
and spd basis versus rc, (b) charge spilling versus the shell of 
neighbors with non-zero overlaps for the spd basis, and (c) for 
the sp basis. 



FIG. 5. (a) Spilling for occupied and first eight bands for 
diamond with a PAOr,. sp versus rc, (b) charge-spilling versus 
the shell of neighbors with non-zero overlaps basis, (c) the 
same for the spilling of the first eight bands. 
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TABLE I. Bases for silicon (sp), diamond (sp), and alu- 
minum (sp and spd). Sq stands for the charge spilling, and S4 
and Ss for the spilling considering the first 4 and 8 bands, re- 
spectively. PAOoo and PAOrc are fixed bases, whereas STOl, 
ST02, ST03, and PAO^' have been optimized according to 
their respective variational parameters, to minimize Sq. Qs, 
Qp, and Qd are the total charges assigned to s, p, and d 
orbitals. The mean square error of the LCAO bands with re- 
spect to the PW reference bands (in eV) is displayed for the 
valence bands (cr^) and two first conduction bands (uc) for sil- 
icon and diamond, and for the first and second bands {ai+2) 
and third and forth ((73+4) for aluminum. Min S corresponds 
to the spilling lower bound (see text). 





Basis 




Sq 


Ss 


Qs 


Or, 

^p 






<7"c 




PAOoc 


sp 


0.0080 


0.1596 


1.32 


2.68 




0.30 


2.15 




PAOS?* 


sp 


0.0078 


0.1621 


1.35 


2.65 




0.29 


1.97 




P AO 5.0 


sp 


0.0139 


0.1591 


1.40 


2.60 




0.36 


2.09 


Si 


STOl 


sp 


0.0109 


0.1764 


1.48 


2.52 




0.38 


2.43 




ST02 


sp 


0.0076 


0.1562 


1.36 


2.64 




0.34 


2.34 




ST03 


sp 


0.0074 


0.1548 


1.34 


2.66 




0.29 


1.95 




Min S 


sp 


0.0044 


0.0689 














Min S 


spd 0.0001 


0.0014 














Basis 




Sq 


Ss 


Qs 


Qp 




(Tv 






PAOoc 


sp 


0.0035 


0.0719 


1.05 


2.95 




0.32 


2.19 




PAOS?* 


sp 


0.0027 


0.0780 


1.07 


2.93 




0.25 


2.61 




PAO4.0 


sp 


0.0041 


0.0672 


1.10 


2.90 




0.23 


2.22 


C 


STOl 


sp 


0.0039 


0.0771 


1.11 


2.89 




0.26 


3.07 




ST02 


sp 


0.0039 


0.0775 


1.12 


2.88 




0.24 


3.05 




ST03 


sp 


0.0024 


0.0790 


1.09 


2.98 




0.17 


3.36 




Min S 


sp 


0.0008 


0.0226 














Basis 




Sq 




Qs 


Qp 


Qd 


^1+2 


0"3-|-4 




PAOoo 


sp 


0.0165 


0.1141 


1.13 


2.87 




2.03 


3.45 




PAOoo 


spd 0.0022 


0.0053 


0.77 


1.39 


0.84 


0.61 


0.24 




PAOS?* 


spd 0.0011 


0.0079 


0.98 


1.41 


0.61 


0.35 


0.84 




PAOe.o 


sp 


0.0157 


0.1012 


1.30 


1.70 




1.54 


2.97 


Al 


PAOe.o 


spd 0.0014 


0.0042 


1.02 


1.39 


0.59 


0.20 


0.45 




PAOb.3 


sp 


0.0158 


0.0930 


1.34 


1.66 




1.34 


2.63 




PAOb.3 


spd 0.0007 


0.0020 


1.15 


1.45 


0.40 


0.97 


2.15 




Min S 


sp 


0.0009 


0.0025 














Min S 


spd 2.9 lO"'^ 


0.0007 













TABLE II. Basis optimizations for different number of 
bands, for Si, C (diamond), BN, and Al. Superindex indi- 
cates the number of bands used for the <S minimization {v 
means the valence bands). The basis is minimal, except for 
Al, for which it is spd. Sq is the charge spilling, and S4, and 
Ss are the spilling considering the first 4 and 8 bands, re- 
spectively. The mean square error of the LCAO bands with 
respect to the PW reference bands (in eV) is displayed for 
the valence bands (a„) and two first conduction bands (uc) 
for Si, C, and BN, and for the first and second bands ((T1+2) 
and third and forth ((73-1-4) for Al. E^^p and Egap are the di- 
rect gap at r and the absolute gap, respectively, both in eV. 
k is the wave vector of the minimum of the conduction band 
in units of n/a, being |k| = 1 for the X point. 

Basis Sq Ss ay ac E^^^ Eggp |k| 



PAOL 
PAOL 


0.0078 
0.0348 


0.1621 
0.1233 


0.28 
0.94 


1.97 
2.13 


2.91 
2.33 


1.69 

1.44 


0.63 
0.80 


PAOg.o 
PAOi.o 


0.0136 
0.0230 


0.1453 
0.1170 


0.44 
0.91 


2.16 
2.47 


2.78 
2.50 


1.96 
0.95 


0.62 
0.83 


STOl" 
STOl® 


0.0109 
0.0207 


0.1764 
0.1541 


0.38 
0.60 


2.43 
2.27 


2.83 
2.40 


2.60 
1.94 


0.47 
0.60 


ST02'' 
ST02® 


0.0076 
0.0306 


0.1562 
0.1179 


0.34 
0.79 


2.34 
1.48 


2.66 
2.39 


2.09 
0.85 


0.53 
0.80 


ST03'' 
ST03® 


0.0074 
0.0465 


0.1548 
0.1121 


0.29 
1.20 


1.95 
1.28 


2.66 
2.52 


1.74 
0.53 


0.60 
0.83 


PW 










2.55 


0.45 


0.85 


Basis 


Sq 


Ss 






^gap 


Egap 


|k| 



c 



Basis 


Sq 


Ss 


(Tv 


(Tc 


^gap 


Egap 


|k| 


PAO^ 


0.0022 


0.0766 


0.24 


3.06 


10.90 


8.56 


0.93 


BN PAOL 


0.0154 


0.0534 


0.96 


2.21 


7.93 


4.63 


0.93 


PW 










8.19 


3.82 


0.93 


Basis 


Sq 


Si 


Cl+2 


f 3-1-4 









PAOSo 0.0027 

PAOi, 0.0192 

STOl" 0.0039 

STOl® 0.0208 

ST02" 0.0039 

ST02® 0.0293 

ST03" 0.0024 

ST03® 0.0335 

PW 



0.0780 0.25 2.61 

0.0600 0.89 2.04 

0.0771 0.26 3.07 

0.0592 0.94 1.95 

0.0775 0.24 3.05 

0.0545 1.37 1.88 

0.0790 0.17 3.36 

0.0526 1.32 1.46 



6.57 5.83 0.50 

5.19 4.13 0.67 

5.78 5.78 0.00 
4.67 4.34 0.42 

5.79 5.79 0.00 
4.49 3.51 0.58 

6.32 6.18 0.42 

4.81 3.96 0.67 

5.28 3.43 0.76 



Al PAOSo 0.0011 0.0079 0.35 0.84 
PAO* 0.0036 0.0024 0.37 0.23 
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TABLE III. Interaction range for minimal sp bases in 
Si. NS stands for the number of neighbor shells with 

non-ncgligiblc interactions and rQg g% is a measure of the size 
of the orbitals as defined in the text. Sq stands for the charge 
spilling, and Sa for the spilling considering the first 8 bands. 
The mean square error of the LCAO bands with respect to 
the PW reference bands (in eV) is displayed for the valence 
bands and two first conduction bands (ctc). STOI1.75 
and STOli.85 are the STOl bases for the exponents P = 1.75 
a.u.~^ (Ref. 5) and 1.85, respectively. 



Basis Sq 


Ss (Ty CTc 


^99.9% 


NS 


STOI1.75 0.0319 


0.2279 0.96 4.27 


5.1 


4 


STOI1.8B 0.0514 


0.2533 1.07 4.64 


4.8 


3 


PAOoo 0.0080 


0.1596 0.30 2.15 


7.9 


9 


PAO5.0 0.0139 


0.1591 0.36 2.09 


4.6 


3 


PAO4.0 0.0506 


0.2269 1.01 3.89 


3.7 


2 


TABLE IV. Singlc-z, double-z (2-z), and s 


'pd for Si, C (di- 


amond). Sq stands for the charge spilfing, and Ss foi 


■ the 


spilling considering the first 8 bands. The mean square error 


of the LCAO bands with respect to the PW reference bands 


(in eV) is given for the valence bands (a-u) and two first 


con- 


duction bands ((Tc). Superscript v means optimized for the 


valence. STO 2-z bases after Ref. 21. 






Basis 


Sq S8(4) 


et, 




ST03'' sp 


0.0074 0.1548 


0.29 


1.95 


STO 2-z 


0.0099 0.0735 


0.36 


1.50 


STOS" spd 

Si 


0.0009 0.0101 


0.09 


0.23 


PAO5.0 sp 


0.0139 0.1591 


0.36 


2.09 


PAO5.0 2-z 


0.0085 0.0827 


0.36 


1.50 


PAO5.0 spd 


0.0042 0.0132 


0.15 


0.38 


ST03" sp 


0.0024 0.0790 


0.17 


3.36 


STO 2-z 


0.0032 0.0059 


0.36 


0.40 


C 

PAO4.0 sp 


0.0041 0.0672 


0.23 


2.22 


PAO4,. 2-;: 


0.0018 0.0105 


0.28 


1.37 



TABLE V. Band gap in Si and C (diamond) for single- 2; 
and double-2 {2-z) bases. Eg^p and Egap stand for the direct 
gap energy at F and for the absolute gap energy, respectively, 
both in eV. k is the wave vector of the minimum of the con- 
duction band in units of tt/o, being |k| = 1 for the X point. 





Basis 




^gap 


Egap 


|k| 




ST03'' 


sp 


2.65 


1.74 


0.60 




STO 


2-z 


2.40 


0.88 


0.73 




ST03'' 


spd 


2.58 


0.76 


0.85 


Si 














PAO^ n 


bp 


2.47 


2.18 


0.47 




PAO5.O 


2-Z 


2.49 


1.03 


0.70 




PAO5.0 


spd 


2.65 


0.91 


0.85 




PW 




2.55 


0.45 


0.85 




ST03'' 


sp 


6.32 


6.18 


0.41 




STO 


2-z 


5.46 


4.30 


0.73 


c 














PAO4.0 


sp 


5.98 


5.89 


0.34 




PAO4.0 


2-z 


6.00 


5.28 


0.53 




PW 




."').2S 




0.76 



TABLE VI. Calculated charge transfer using optimized 

PAO^g* basis for some zincblende semiconductors. Qc and 
Qa stand for the valence charge on the cation and the anion, 
respectively, 5Q for the charge transfer with respect to neutral 
atoms, and Sq for the charge spilling. Numbers in parenthe- 
sis were obtained from Hartree-Fock LCAO calculations,^^ 
except for OiuVs. for wliicli LD.V LCAO was used." ' 





Basis 


Sq 


Qc 


Qa 


SQ 


BN 


sp 


0.0022 


2.19 (2.14) 


5.81 (5.86) 


0.81 (0.86) 


BP 


sp 


0.0038 


3.51 (3.34) 


4.49 (4.66) 


-.51 (-.34) 


AlP 


sp 


0.0035 


2.15 (2.20) 


5.85 (5.80) 


0.85 (0.80) 


SiC 


sp 


0.0071 


2.30 (2.19) 


5.70 (5.81) 


1.70 (1.81) 


GaAs 


sp 


0.0041 


2.58 


5.42 


0.42 




spd 


0.0010 


2.78 (2.88) 


5.22 (5.12) 


0.22 (0.12) 



TABLE Vll. Calculated excess charge, in units of e, for 
the upper and lower Si atom in the dimer of a Si(100)-2xl 
surface. A minimal sp basis has been used in all cases. Sq 
stands for the charge spilling. PAO^?* and STOl have been 
optimized for bulk silicon. STOI1.75 is a STOl basis with 
13 = 1.75 a.u.-^ (Ref. 5). 



Basis Sq SQup &Qdo-wn 

PAOoo OO8O ^OlO -f-0.12 

PAOS?* 0.0078 -0.10 -F0.13 

PAOs.o 0.0139 -0.13 -hO.lO 

STOli.TB 0.0318 -0.15 -hO.lO 

STOl 0.0076 -0.13 -1-0.12 
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